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I. INTRODUCTION 

Recent experimental progress in the field of ultracold 
atomic gases loaded in optical lattices 1 opened new per- 
spectives in the research of isolated quantum systems 
out of equilibrium 2-6 . In particular it has made possi- 
ble to study the evolution of many-body strongly cor- 
related closed systems after a quantum quench, i. e., 
an abrupt change of one of the parameters of the sys- 
tem's Hamiltonian. Such out-of-equilibrium setups can 
be cleanly realized in the experiments with cold atoms 
due to the fine control over the effective parameters of 
the system, such as external fields and interactions be- 
tween the atoms, even in real time. In parallel, there has 
been an intense theoretical effort in the study of the out of 
equilibrium dynamics in many-body systems (see Ref. 7 
for a review). The investigations have focused mostly on 
the features of the steady state after a quantum quench, 
the conditions under which it is described by a (gener- 
alized) Gibbs ensemble and the existence of metastable 
states at intermediate times and whether they relax for 
longer times. Other relevant and interesting aspects that 
arise are related to the dynamics after the quench in con- 
nection to the so called 'light-cone effect' 8,9 . 

These aspects of the non-equilibirum dynamics after 
a quantum quench have been investigated in various 
specific models such as the quantum Ising chain 10-14 , 
one dimensional (ID) bosonic models 15-18 , the sine- 
Gordon 19-21 , and the Luttinger model (LM) 22 < 23 . Gen- 
eral results have been obtained from theoretical investiga- 
tions involving conformal field theory (CFT) 9,24 ' 25 . The 
LM, and other closely related ID models, are among the 
most attractive candidates for studying ID systems out 
of equilibrium due to the availability of a manageable 
exact solution from which one can extract correlation 
functions in closed form. On the other hand, the (short- 
range) LM is the fixed point of a large class of ID gapless 
systems, called Luttinger liquids; the question of under 
what conditions this universality property is extensive to 
the non-equilibrium domain constitutes an open problem. 



A priori we cannot expect the results obtained for the 
LM out of equilibrium to hold for more general systems, 
since couplings that are irrelevant in equilibrium, might 
become relevant when the system is driven out of equi- 
librium 26 , at least for long times. From another point of 
view, the dynamics after a quenching process will excite 
high-energy states that are not well described by the LM. 
Still, in Refs. 27-29 it has been shown that some char- 
acteristics of the long-time dynamics after a short-range 
interaction quench of certain ID systems are universal in 
the Luttinger liquid sense. 

So far, theoretical investigations on quench dynamics 
have focused primarily on systems that interact via con- 
tact or, more generally, short-range potentials, i.e., that 
decay faster than the inverse of the distance between par- 
ticles. The origin of this bias can be found in the short- 
range character of the couplings in the prevailing experi- 
mental setups that motivated the field in recent years 2,3 . 
Even in the more traditional condensed matter systems, 
the Coulomb interaction between electrons is very often 
screened, particularly in non-isolated ID systems 30 , giv- 
ing rise to an effective short-range interaction description 
in most of the cases. However, this does not exhaust the 
possibilities of studying out of equilibrium systems, and 
in particular the consequences of a long-range potential 
seems to be unexplored, though some results are available 
in a Hubbard chain with long-range hopping 31 . A con- 
crete realization of particles interacting via a long-range 
potential are linear ion traps 32 . These systems, as their 
neutral atomic counterparts, turn out to be highly tune- 
able and are able to maintain coherence for long times, 
becoming thus ideal candidates for studying out of equi- 
librium phenomena. 

Regarding the models used in the past to account 
for long-range potentials, it is specially interesting the 
case of the LM with Coulomb interactions. Its equilib- 
rium properties have been studied in a series of theo- 
retical works 33-35 , leading to some well known results, 
such as the prediction of the Wigner crystal phase for- 
mation of ID electrons, which has been recently mea- 
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sured in isolated carbon nanotubes 36 . From a perturba- 
tive perspective, the fast exponential decay of the (fourier 
transformed) long-range potential with the momentum 
exchanged between colliding particles ensures that the 
sudden connection of interactions is likely not to scatter 
particles to high-energy states. Therefore, we expect the 
Coulomb LM to faithfully capture universal properties of 
more general long-range systems. 

In this work we shall consider the non-equilibrium dy- 
namics after a quantum quench of the LM with long- 
range interactions. We pay most of the attention to the 
case of a sudden turning on of the Coulomb repulsion 
between the fermions starting from an arbitrary short- 
range interacting initial state, though we briefly discuss 
the opposite case, of a sudden turning off of the Coulomb 
potential. We analyze in detail the finite time dynam- 
ics after the quench, and show, in particular, that cor- 
relations spread with increasing velocity, allowing for a 
propagation of information faster than in systems with 
short-range interactions. 

The rest of the paper is organized as follows: in section 
II, we introduce the LM and describe its exact solution in 
terms of bosonic quasiparticles (QPs). The main results 
of the work are presented in section III, where we de- 
scribe the stationary state correlations and the finite time 
dynamics. We consider the case where the long-range in- 
teraction between the fermions is suddenly switched on 
and the reverse situation, that is, when the interaction 
in suddenly switched off. In section IV we proposed an 
experimental realization of the studied system using ul- 
tracod ions loaded in electromagenetic cavities. Finally, 
we summarize the main results in section V. 



II. THE LUTTINGER MODEL AND THE 
BOSONIZATION SOLUTION 

The LM describes a system of interacting fermions in 
one dimension 37-40 . The key simplification of the model 
is the assumption of a linear dispersion relation for the 
free fermions, which induces a clear separation between 
right and left moving particles. The LM Hamiltonian is 

Hlm = H + H 2 + H4, (1) 

where 

H = ^2 v F q : c\ r c q , r :, (2) 

q,r=R,L 

is the free piece of Hlm- Here c q r and c q:T are fermionic 
creation and annihilation operators at momentum q and 
vf is the Fermi velocity. The index r labels the chirallity 
of the particles, denoting left (L) and right (R) moving 
fermions. The fermionic normal ordering denoted by : 
... : is needed to remove from the expectation values 
the infinite contribution arising from the fact that the 
ground state is a Dirac sea 41 , namely, a state where all 
single-particle fermion levels with q < are occupied. 



The strength of the scattering can be parametrized us- 
ing two functions, 172(9) and 54(g), that are related to pro- 
cesses that exchange fermions between the two branches 
of the spectrum and those that leave the fermions on its 
original branch respectively. In terms of these, the inter- 
acting piece of Hlm reads 

H-2 = Y^g2{q) ■■ PqRpqL (3) 
HA = ~ ^ 9i(q) ■ PqRP-qR + PqLP-qL ■, (4) 

<? 

where we have defined densities in momentum space as 
Pgr = Sfc • c \-qr Ckr '•• ^ n wna t follows we will take 
52(9) = .94(9) = V(q) as the Fourier transform of the two 
body interaction potential between the fermions V(x). 

The elementary excitations of the system are bosonic 
QPs describing low energy collective density modes of the 
system 39,40 . In order to see this, we first note that the 
density operators p qr obey the following commutation 
rules: 

[p\riPr'q'\ = -rS rr > 6 qq , 7l q , (5) 

where we have defined the integer n q = The spe- 
cial algebra of Eq. (5) can be transformed into the usual 
bosonic commutation relations by introducing the follow- 
ing operators: 

b{q) = _L= ( '«*' *» q I ?' (6) 

Notice that the 9 = components (the zero modes) re- 
quire a separate treatment since po r — N r is the de- 
viation, relative to the ground state, in the number of 
fermions of a given chirality. It is customary to intro- 
duce the combinations 

N = N R + N L , J = N r -N l . (7) 

In therm of these, and of the bosonic operators intro- 
duced in Eq. (6) the Hamiltonian of the LM writes 

H = J2v F \q\blb q + ^(N 2 + J 2 ), (8) 

q^O 

H * = lj2 V ^)\Q\hb-q + blbl q ] (9) 

+^V 2 -J 2 ) (10) 
H 4 = *£V(q)\q\blb q (11) 

q^O 

+ ^V 2 + J 2 ). (12) 

Ignoring the zero mode sector (that will not contribute 
in the L — > oo limit), we can diagonalize the above 
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bosonic Hamiltonian by means of a Bogoliubov trans- 
formation by introducing new operators a q and a q as fol- 
lows: 



,t 



cosh #((7) s'mhd(q) 



aL q J \ smhO(q) cosh8(q) J y b l _ q 



(13) 



where the parameter of the transformation 6(q) satisfies 
the relation 



tanh 26»(g) 



V(q) 



2nv F + V(q)' 



(14) 



In terms of the new operators the Hamiltonian is ren- 
dered diagonal: 

Hlm = Hq + H2+ H4 = e(q)a q a q +zero mode terms, 

q^Q 

(15) 

where e(q) — Vp\q\\Jl + 2-^1 i s the dispersion relation 

of the bosonic QPs. This defines the equilibrium solution 
of the LM, with which all correlation functions can be 
calculated. This is of particular importance since the 
equilibrium LM is the renormalization group fixed point 
of a large class of gapless one dimensional models 41 ' 42 , 
called Luttinger liquids. 

For systems with short-range interactions the LM cal- 
culations give rise to the characteristic power law de- 
cays of the correlations 40 ' 41 . It is worth clarifying that 
in this work we shall call short-range potential to any 
potential V(x) whose Fourier transform V(q) is finite at 
q = 0, falling in this category potentials that are usu- 
ally termed as long-ranged, such as the dipole-dipole in- 
teraction V(x) ~ a; -3 . In such cases, the asymptotic 
form of several correlation functions is controlled by the 
value of the potential at q = 0, through the so-called 
Luttinger parameter K = (1 + ^7) ^ 2 and the renor- 

malized velocity v = Vp yl + ■ For systems in which 

the Coulomb potential is unscreened (V(q) diverging log- 
arithmically for q — > 0, see Section III A) the LM predicts 
the existence of the analogous of a Wigner crystal (the 
crystalline phase of electrons) in ID. In particular, as was 
shown by Schulz 33 , the Ak F part of the density-density 
correlation function decays as e~ clogl/2 ^ x \ much slower 
than a power law, signaling the stability of a phase in 
which the electrons fluctuate around equidistant lattice 
sites. 



III. INTERACTION QUENCHES IN THE 
LUTTINGER MODEL 

Even though general quenches between quadratic 
Hamiltonians such as the LM can be fully solved 22 > 23 ; 
we shall focus on a specific situation in which the range 
of the interaction is suddenly modified. Let us assume 
that at t = the system is prepared in the ground state 
\GS)i of the Hamiltonian Hi, defined through Eq. (1) 



with an interaction potential Vi(q). The Hamiltonian 
that henceforth dictates the temporal evolution, Hf, has 
a different interaction potential Vf(q). At the level of the 
bosonic representation, both Hamiltonians can be put in 
diagonal form by means of canonical Bogoliubov transfor- 
mations characterized by different parameters 9i(q) and 
0f(q), related to the respective potentials by Eq. (14). 
The solution of the problem is then reduced to a change 
of basis consisting of a series of nested Bogoliubov trans- 
formations performed in order to find the temporal de- 
pendence of the operators that diagonalize the bosonic 
version of Hamiltonian Hf. We show the details of this 
procedure in the Appendix A. We shall mostly consider a 
simple protocol, according to which the interactions de- 
scribed by H 2 and H4 are suddenly changed from short 
to long-range, being the case of an initial non-interacting 
Hamiltonian (32(9) = gi{q) = 0) a specific example. We 
shall deliver most of the attention to this protocol, and 
briefly analyze the reverse one in Sec. IIIC. 



A. Turning on Coulomb interactions 

Consider that initially we have a one-dimensional gas 
of electrons interacting via a short-range potential, whose 
subsequent temporal evolution occurs with the electrons 
interacting via the Coulomb potential. Since the inclu- 
sion of a finite but short range initial potential would 
only renormalize the value of the Luttinger parameter 
Ki, we shall only consider a contact potential V%{x) — 
s5{x), which induces an initial Luttinger parameter Ki = 
(1 + ^r^r) -1 ^ 2 - We can always recover the limit of a 
non- interacting initial state by setting Ki = 1. For the 
Coulomb potential we take the form Vf(x) = ^ x 2 +d2 

with Fourier transform Vf(q) — 2e 2 K (qd), where -Ko(C) 
is the modified Bessel function of the second kind of order 
zero and d is a short-distance cutoff. This phenomeno- 
logical form has been proposed to effectively describe the 
unscreened Coulomb potential on semiconductor wires 43 , 
in which case d is associated to the transverse dimensions 
of the wire. The spectrum of the bosonic QPs of the 
Hamiltonian is of palsmonic type in ID: 



= \q\vp^l + 2gK (qd) ~ \q\ log 



1/2 



(16) 



with g = , and the logarithmic form is the approxi- 
mate result for small q. The energy of the QPs goes to 
for q — > 0, nevertheless its group velocity Vf(q) = d<Li dq 

diverges as vpyjl — 2glog(qd) for q <C d^ 1 . 

In order to gain insight into the properties of the sys- 
tem following such a quench, we will compute static 
(equal time) correlation functions of the form Cq(x, t) = 
(e lHft O(x)O(0)e~ lHft ), where O is an operator and (. . .) 
stands for the average taken over the initial state. Notice 
that, since in general the initial state is not an eigenstate 
of Hf, time translation invariance is broken and Cq{x, t) 
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is explicitly time dependent. To begin with, let us con- 
sider the one-particle density-matrix 



-iHft\ 



(17) 



whose Fourier transform leads to the instantaneous mo- 
mentum distribution function. This correlation function 
can be calculated using the bosonization identity 



\/2na 



itp r (x) 



(18) 



where 

(f> r (x) 



q=£0 



-aq/2 



(e ir ^b qr + e- ir " x bl r ), (19) 



a is an UV cutoff used to regularize the short distance 
divergences of the model and the F r , usually called Klein 
factors, are unitary operators obeying {F£,F r '} = 25 rr i . 
After the quench we obtain the following result for this 
correlation function in the thermodynamic limit (see de- 
tails in Appendix B): 



C^{x,t) = cf r {x) 



(Ki-K- 1 ) 



$ 2 (x,t) 



(20) 



where the functions 



dq 



e~ aq sinh 2 (29 /(<?)) 



x sin 2 (e f (q)t) (I - cos(qx)), (21) 



and 



9 2 (x, f ) = 2 J o y e ~ aq cosh(29 f (q)) 

x sinh(20/(g)) sm 2 (e f (q)t)(l - cos(gx)) (22) 

are independent of the initial interactions. In these equa- 
tions 



v Z7T 



(23) 



is the correlation function at t = 0, i.e., the equilibrium 
equal-times correlation function for a system described 
by Hi. 



1. Steady state and the generalized Gibbs ensemble 

In order to investigate the steady state properties af- 
ter a quench to a Coulomb potential, we note that for 
t — > oo the oscillatory, time-dependent factor within the 
integrals in Eqs. (21) and (22) average to 1/2, regardless 
of the divergence in the mode velocity. By evaluating 



the resulting integral we obtain the long distance i>(! 
behavior of the one-particle density matrix, which takes 
the following asymptotic form: 



C^ r (x,oo) ~ C§l(x)e 



\Ki \og 2 {x/d) 



(24) 



We thus see that the leading order is given by the correc- 
tion introduced by the quench, that decays faster than 
any power law, and in particular decays faster than the 
equilibrium correlations in the system with long-range 
interactions, where it has a leading behavior given by 
e -eiog 3/ (x/d)33_ j n ac yition to the presence of C^\x) as 
an overall prefactor, the memory of the initial condition 
is also reflected by the factor that enters in the ex- 
ponent of the correction introduced by the quench and 
modulates the coupling g. 

Few-body correlations after a quantum quench in many 
exactly solvable models, in particular in those that 
can be factorized in separate momentum sectors, have 
been shown to relax to a generalised Gibbs ensemble 
(GGE) 16 ' 25 ' 44,45 . The LM belongs to this group, even 
though some other quantities such as energy fluctuations 
cannot be obtained from the GGE 22 ' 23 . In Ref. 46 two 
main ingredients were identified as necessary to show the 
equivalence of the steady state static correlations with 
the GGE, namely, dephasing among different Fourier 
components and noncrgodicity of the correlations, i.e. 
the fact that asymptotically correlations depend only 
on the eigenmode occupations. Other integrable mod- 
els have been shown to relax to the GGE 28,47 , though 
performing the trace associated to the GGE for specific 
interacting systems represents an almost impossible task. 
Moreover, a general result, which is always desirable, re- 
lating relaxation in integrable models to the GGE is still 
lacking. 

In a quantum quench into the Coulomb Luttinger 
model, the singularity present in v(q) when q — > does 
not affect the dephasing (since the energy remains finite 
for any q) while preserves the nonergodicity of the corre- 
lations implying therefore the equivalence with the GGE. 
However, the equivalence of dynamic correlations is less 
clear. Let us show that the equivalence holds, by consid- 
ering as an example the stationary limit of the two times 
Green's function: 

C^ r (x,t,ta) = ($(*»* + *oM-(0,to)), (25) 

where, as before, (. . .) stands for the expectation value 
taken over the initial state, that we will take as a non- 
interacting Fermi gas (Ki = 1) for simplicity. After us- 
ing the bosonization identity Eq. (18) and the Baker- 
Hausdorff formula, we obtain: 



C^(x,t,t ) - 1 !_ e lM*,t+to),Mo,to)) 



L 



x e 



i{<j> r (x,t+to)-<f>r-(0,t o ))\ 



(26) 



where the commutator is a c-number that only depends 
on t. Since we are dealing with a Gaussian theory, we 
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then make use of the property (e J 



-5<^ 2 ). The 



quantity of interest is essentially the exponent, 
((<M:c,< + to)-0r(O,io)) 2 } = 



q>0 



-{\e lqX f(q,t + t )-f(q,t Q )\ 2 
+ \e-^g(q,t + to)-g(q,t )\ 2 } (27) 



which, by using the expressions for f(q,t) and g(q,t) 
given in Appendix A can be shown to reduce in the limit 

to — > oo to 



in a system with short-range interactions. The fact that 
in our system the dynamic correlators in the steady state 
can be mimicked by a GGE suggests that the short-range 
interaction hypothesis can be relaxed in some cases. This 
is related to the fact that although there is no finite 
speed limit for the information propagation carried by 
the plasmonic QPs, correlations change their behavior 
when crossing a non-linear light cone, as we will discuss 
in the next section. 



lim ((4> r {x,t + t )-4> r (0,t )) 2 ) = 

Zq—>OQ 



q>0 



- aq 



x {cosh 2 (29 f(q)) (1 — cos(qx) cos(e/(g)f)) 
cosh(20/(g)) sin(gx) sin(e/(g)t)}. (28) 



Next, we introduce the GGE with which we shall com- 
pare the above result. In terms of the bosonic basis that 
diagonalizes the evolution Hamiltonian Hf, {a q , a q }, (see 
Appendix A), the first natural choice for the GGE den- 
sity operator is 



1 

PgG = -w— exp 



q>0 



(29) 
, and 



where n(q) = a\a q , Z gG = Tr exp J2 q> o Hq)n{q) 

the Lagrange multipliers X(q) are fixed by the initial con- 
ditions: 



A(<?) = (n(q)) eG = (n(q)) = sinh 2 ^ (q)), 



(30) 



where (0) g c = Tr[p g QO]. In the following we will argue 
that, if we define Cf^(x,t) = (ipfa, t)i/) r (0, 0)) gG , then 

lim C ilr (x,t,t )=Cf(x,t). (31) 

to— >co ^ r 

To show this, we proceed in complete analogy with the 
previous calculation. We obtain in the GGE the result 



Cf r {x,t) 



(1 



-2ira/L\ 



,[0,(x,t),<MO,O)] 



L 



x exp 



; ((<Mx,t)-<Mo,o)) 2 ) gG 



(32) 



Again we are only interested in the exponent, being the 
remaining factors equal to Eq. (26). A calculation anal- 
ogous to the one leading to Eq. (28) shows a complete 
equivalence with the exponent of Eq. (32) . 

Interestingly, it was recently shown 48 that the same 
GGE that reproduces the static (one time) correlations 
after relaxation can describe the dynamic (two times) 
correlation functions. The proof is based on the Lieb- 
Robinson bounds 8 on the information propagation speed 



2. Dynamics and non-linear light cone effect 



Now we turn to the subject of finite time dynamics of 
the system after the quench. Generally, after a quan- 
tum quench the space-time dependence of the correla- 
tions in systems with short-range interactions exhibit the 
so-called light-cone effect: in the region of space such that 
x > 2vt, where v is the characteristic velocity of the ex- 
citations, correlations essentially keep the form already 
present in the initial state. In the complementary region 
t > x/2v, they are of the form of the steady state cor- 
relations. This effect has been rigourously demonstrated 
for one dimensional models with conformal invariance 9 
of which the LM is a specific example, and can be un- 
derstood as follows: in these models the dispersion rela- 
tion of the QPs is linear e(q) — vq, with v = the 
(constant) group velocity of the QPs. The initial state (a 
complicated excited state of the Hamiltonian Hf) acts as 
a source of QPs that propagate semiclassically through 
the system with the velocity v. A pair of QPs emerging 
from the same point is entangled, and their arrival to two 
distant points affects the correlations initially present be- 
tween them. This picture has been validated in a series of 
theoretical works 19,22 ' 23,49-52 and even experimentally 4 . 
This simple interpretations faces a problem in one di- 
mension when the interactions are long-ranged, since the 
group velocity Vf(q) diverges for q — > whereas the en- 
ergy is finite [see Fig. (1)] . Moreover, the QP momentum 
distribution (a q a q ) = sinh 2 (6 f (q)) (using the notation of 
the Appendix A), which is a constant of motion, is peaked 
around q — 0, the region where the velocity varies more 
abruptly. Thus, if we stick to the above described pic- 
ture, we would expect that the correlations of the steady 
state propagate instantaneously throughout the system. 
We shall see that this is not the case and point out the 
failure of the above reasoning when long-range interac- 
tions are present. 



We first show some analytical results in the asymp- 
totic regimes. We start by analyzing the integral defin- 
ing $i(x,t) in Eq. (21) (being the analysis completely 
analogous for $2(^1 1)) since this integral define the cor- 
rection induced by the quench on top of the initial condi- 
tion C^\x). After integration by parts we find that for 
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FIG. 1: (Color Online) Spectrum e(q) (solid line) and group 
velocity v(q) (dashed line) of the bosonic QPs of Hf, for g = 1. 
e(q) is given in units of ^f- and v(q) in units of vf- 




x d and fixed t 



where 



2 log(dx _ 1 ) sinh 2 (20/ (x~ 1 ) ) sin 2 (e/ (a: - 1 )t) 

+ t(t), (33) 



£(t) - 2 f ^ sinh 2 (20/(9)) sin 2 ( e/ («)f) 

~-|bg 2 (2« J rt/d) (34) 

is a function that only depends on t. The last equality is 
valid for vpi ^> d 68 . The first term in Eq. (33) vanishes 

approximately as — ^ for 2ef(x~ 1 )t -C 1, represent- 
ing a subdominant correction with respect to the corre- 
lations in the initial state, Eq. (23). The same analysis 
holds for $2(2;, t). We thus see that the correction in- 
duced by the quench on the one-particle density matrix 
reduces, up to exponentially suppressed corrections, to a 
time-dependent factor: 



C^t) = cf(x)e-i K ^ 2 ^ t ' d \ 



for times such that 
d 



< t < t x = . 

vf 2v F \/l - 2glog(da 



(35) 



(36) 



Here we have used the explicit form of the plasmon dis- 
persion Eq. (16), and we consider distances x 3> d. For 
such short times, spatial correlations have the same form 
as in the initial state, i.e. are given by C$(x). 

Note that the denominator of Eq. (36) is essentially 
the group velocity 2vf(x~ 1 ) and therefore we can rewrite 
Eq. (36) as 2v r^=r\ ^ t. For a short-range potential 
the group velocity Vf(q) tends to a constant v/(0) as 
q — >• 0, thus determining a linear light-cone 2 v(fi) = * 
for sufficiently large x. On the contrary, the slow group 
velocity divergence at q — > of the long-range potential 



generates a (weak) non- linear light- "cone" as suggested 
by Eq. (36), in which the behavior of the correlation 
function crosses over from a short times regime to a long 
times one. This result can be understood in terms of the 
following picture: the peculiarities of the ID plasmon 
dispersion reflect in that for decreasing q, the energy of 
the plasmonic modes decreases whereas its velocity di- 
verges [see Fig. (1)]. Consequently, and since a mode 
with momentum q can only propagate correlations (infor- 
mation) over distances larger that 1/q (its wavelength), 
the minimum distance over which the new correlations 
can propagate is greater for the faster modes. Alterna- 
tively, at short times only high energy modes are active, 
but these have low velocity and therefore do no propa- 
gate far. At larger times, the lower energy modes that 
get progressively activated are faster and therefore reach 
more distant points. This argument is in agreement with 
Eq. (36). Hence, long-range interactions reveal more ex- 
plicitly the quantum character of the QPs propagation. 
The sublinear behaviour of the light-cone also implies 
that quantum information can propagate faster than any 
finite speed of sound. However, this effect is not due to 
an acceleration of excitations with defined momentum, 
since each mode has constant velocity v(q), but to the 
quantum nature of the collective behavior, which only 
activate fast modes at later times. Whether this effect 
will be also present in actual lattice models with long- 
range interactions remains an open problem. The fact 
that it lyes on the peculiarities of the plasmon dispersion 
in ID suggest its generality. In such a case, and since 
deviations from linearity in the light cone are not strong 
but just a logarithmic correction, simulations with cur- 
rent algorithms should be feasible. 

Cases of non-linear information propagation have also 
been reported in systems with short-range interactions, 
however, the mechanisms behind are totally different. 
For instance there has been detected certain type of 
bosonic models for which the light cone can "bent out- 
wards" (accelerated information propagation) 53 . In addi- 
tion, it has been predicted that in disordered spin chains 
the light cone can "bent inwards", with the radius of the 
light cone growing logarithmically with time 54 ' 55 . 

In order to complement the previous discussion we 
shall present some numerical calculations. First we will 
carefully analyze the behavior of &i(x,t) in the case 
where we suddenly turn on short-range interactions. In 
Fig. (2) we show the temporal dependence of $1(2:, t) for 
different values of x after turning on a contact potential 



V(q) = V, 
and a gaussian potential 

V(q) = V(0) exp 



g 

2f7 2 



(37) 



(38) 



At short times, initial state correlations dominate and 
correlations follow the x — > 00 curve. Then deviate from 
the asymptotic behaviour, develope a maximum and de- 
cay to the long-times value signaling the stationary state 
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FIG. 2: (Color Online) Time evolution of &i(x,t) from nu- 
merical integration for various values of x after suddenly turn- 
ing on (a) a contact potential (Kf — 0.5) and (b) a gaussian 
potential (2a 2 = d~ 2 ). The black dashed curve is, in both 
graphics, the result for x — > oo. The vertical lines are, for 
each x, the times defined by Eqs. (39) and (40) respectively. 



behavior. The differences between these two cases are 
also interesting. For the contact potential the maximum 
is a sharp peak (that translates into a singularity in the 
derivative of the correlation function) located exactly at 



t = 



2t>' 



(39) 



with v — t>Fy 1 + the group velocity of the collec- 
tive modes. For the gaussian potential the peak broadens 
and oscillations appear mounted on the final decay to the 
constant. The position of the peaks is approximately lo- 
cated at 



t 



2«(0)' 



(40) 



where v(q) is the group velocity induced by the gaus- 
sian potential, but deviations from this value are visi- 
ble. This illustrates our previous argument, according to 
which a finite range potential is related to a linear light 
cone in virtue of the finiteness of the potential at q = 0. 
Moreover, we can conclude that a finite dispersion in the 
velocity distribution of the QPs that carry the new cor- 
relations induces oscillations of the correlation functions 
and broadens the crossover region in which correlations 
change their behavior. 
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FIG. 3: (Color Online) Time evolution of &i(x, t) for vari- 
ous values of x after a quench from a short-range interacting 
initial state to a long-range interacting one with g = 1, from 
numerical integration. The dashed line represents the func- 
tion £(t) defined in Eq. (34). 
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FIG. 4: (Color Online) Density plot of the function $i(x,t) 
for g = 1 from numerical integration. The black dashed line 
is the curve t = 1.26t x (Eq. 36). The solid red line marks the 
positions of the maximums. Region I is the region where the 
initial state correlations dominate. In region II correlations 
have the form of the steady state. 



For the Coulombian potential, Fig. (3), some simi- 
lar qualitative features are found, although in this case 
the peak is appreciably more broader than in the case of 
the gaussian potential since the velocity of the plasmons 
ranges over an infinite set of values. Nevertheless a much 
more interesting difference arises since the maximums do 
not lie on a straight line in the space-time plane (as was 
the case in the previous examples). This can be seen in 
Fig. (4), where we summarize the non-linear light cone 
effect picture. In region / the correlations are the same 



as in the initial state up to a time dependent prefactor 
(Eq. (35)) while in region II the steady state correla- 
tions Eq. (24) dominate. In between / and II there is a 
broad crossover region where the function reaches a max- 
imum and oscillates. The center of the crossover region 
is well estimated by the position of the maximums, that 
is accurately described by the curve t = 1.26t x (Eq. 36). 
Deviations from linearity are weak since they are given 
by a factor log 1 / 2 x/d. 

Finally we note that the QP dynamics has shown to 
be intimately related to the growth of the entanglement 
entropy. For example, in a CFT, the well defined prop- 
agation velocity of the QPs is directly related to the lin- 
ear growth of the entanglement entropy in time after a 
global quench 24 . Therefore, we can speculate that the 
entanglement entropy growth after a global quench in 
the Coulomb LM will follow the same non-linear grow as 
found earlier for the light-cone. However, this can be only 
confirmed through an explicit calculation of this quantity 
which is out of the scope of the present article. 



3. Dynamics after a quantum quench from a 
non-interacting initial state 

Let us analyze next in more detail the case of a non in- 
teracting initial state, for which Eqs. (20)-(21) hold pro- 
vided that one takes Ki — 1. We shall consider two 
quantities that were utilized in studying the universal- 
ity of the Luttinger liquid description after an interac- 
tion quench 27 , namely, the time evolution of the Z-factor 
jump of the instantaneous momentum distribution func- 
tion n(k,t) at the Fermi momentum kp and the kinetic 
energy per length eki n (*)- The discontinuity of the mo- 
mentum distribution is 



Z(t) = lim n(k,t) — lim n(k,t), 



(41) 



where n(k, t) is the instantaneous Fourier transform of 
the one-particle density matrix Eq (17). Z(t) can be 
interpreted as a time-dependent 'Landau quasiparticle" 
weight in an effective Fermi liquid description of the 
system at finite times 22,23 . For quenches involving a 
Hamiltonian with only finite range interactions, and 
thus characterized by parameters Kf and Vf, Z(t) ~ 
t~^ K ~f~ K f ~ 2 ) 23 >56 The same computation in the case 
of a long-range Coulomb potential yields the result 



Z(t) 



■ \og 2 (2v F t/d) 



(42) 



which decays faster than any power law but still slower 
than an exponential law. Thus, at least at the level of the 
one-particle density matrix and the momentum distribu- 
tion, correlations tends to the steady state value faster 
for a quench involving long-range interactions than for 
short-range ones. The behavior of Z(t) is depicted in 
Fig. (5). 
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FIG. 5: (Color Online) The Landau quasiparticle weight Z(t) 
for different values of the coupling constant g, from numerical 
integration. Note the log-log scale; straight lines would be 
expected for power laws. 



We shall also focus on another quantity, the kinetic 
energy per unit length, defined as 

iHft tj „—iHft\ 



2tt 



dq qsmh 2 (29 f {q))sm 2 {e f (q)t). (43) 



It describes the rate at which the excitations of the ini- 
tial state fade away. In the case of turning on finite range 



interactions 



| de kin (t) | 



\{K 2 -KJ 2 -2)v F 



-327 



in the long- 
times limit. For of long-range interactions, the contri- 
bution coming from the divergence of the potential only 
represents a logarithmic correction to the power law, and 
the same t~ 3 law holds, as shown in Fig. (6). We thus see 
that the long-range character of the interactions affects 
some of the aspects of the dynamics, while contributes 
as a subdominant correction to others. Strong evidence 
has been provided 27 pointing to the universality in the 
Luttinger liquid sense of the asymptotic time dependence 
of these two quantities for short-range interactions. The 
extent to which this universality is maintained for long- 
range interactions constitutes an open problem. 



B. Vertex operators correlations and 
particle-number fluctuations 

Other correlation functions exhibit a similar behavior 
to the one-particle density matrix. Let us now focus on 
the correlators: 

jHm[4>(x,t)-tl>(p,t)]\. 



where 



(j>{x,t) 



= e -2m 2 ([0(x,t)-0(O,t)] 2 )^ 
_ /An[e(x,t)-8(0,t)]\ 



e 2 



([$(x,t)-e(o,t)] 2 ) 



<f>R(x,t) - (f)L{x,t) 



(44) 
(45) 
(46) 

(47) 
(48) 
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FIG. 6: (Color Online) The derivative of the kinetic energy 
per unit length for different values of the coupling constant g, 
from numerical integration. In the inset we can see that the 
kinetic energy per unit length reaches a maximum and goes 
to a stationary state constant value, ejdn(i) is given in units 
of 



and 



9{x,t) 



<t>R(x,t) + <f>i,(x,t) 



(49) 



is the canonically conjugated field of <f>(x,t), i.e. 

[<t>{x, t), -d x ,6{x', t)] = i8{x - x'). (50) 

7T 

These correlation functions are related to the total den- 
sity and current-density fluctuations in the LM and in 
Luttinger liquids in general. In fact, in the LM we can 
write the fermionic particle density in real space as : 



d x <j>{x,t) = -n[p R (x,t) + p L (x,t)], 



(51) 



in the L — > oo limit. 

More generally, in the context of the hydrodynamic 
approach of ID interacting systems (the Luttinger liquid 
theory 42 ), C™(x) and Cg{x) describe the wave number 
and phase fluctuations near 2mkp and nkp, respectively, 
of the density correlation function (p(x)p(0)). 

In the Coulomb Luttinger liquid in equilibrium the 
slow decay of the 4k p component of the density 
C?(x) e~ clog 1 x , slower than the 2kp component, hall- 
marks the appearance of the Wigner crystal phase. After 
a quantum quench from an initial state with short-range 
interactions, we obtain in the steady state the following 
correlations: 



Cf(x,oc) ~ cf m {x)(x K ^)' 



g 2 (x/d) 



72 



(52) 
, (53) 



where Cf' m (x) = \x\- 2m * K > and cf ,n (x) = \x\~^ K " 
represent the initial condition. The behaviour of the 
Co correlator is similar to the one encountered for the 



Green's function. Hence, at long distances the correla- 
tions are dominated by the correction factor introduced 
by the quench, that decays faster than any power law, 
and in particular faster than in equilibrium. The initial 
condition is only reflected in the presence of the factor 
Ki multiplying the coupling constant g. On the other 
hand, the correlation function C™ presents more inter- 
esting features: the effect of the quench is, to leading 
order, only reflected in a correction of the initial condi- 
tion's power-law exponent, leading to a decay of the form 
C^(x, oo) ~ \x\~ m Ki for sufficiently large distances. We 
notice that the consequences on C™ of suddenly changing 
the range of the interaction are similar to the ones orig- 
inated in an abrupt change on the value of the stiffness 
parameter K. 

The quantity ([(j>(x, t) — 0(0, t)] 2 ) is also directly related 
to the particle number fluctuations F{x, i) in a segment 
of size x. In fact, by integrating Eq. (51), we obtain that 

F(x,t) = ([N x - (Ay] 2 ) = 7r 2 ([0(x,t) - 0(0, t)] 2 ), (54) 

where N x is the operator that represents the number of 
particles in such region, and T(x, t) is the quadratic mean 
deviation of that particle number. Notice that C™(x,t) 
is the generating function of the cumulants of the (gaus- 
sian) particle number probability distribution. 

The fluctuations in the number of particles of a sub- 
system are intimately related to its entanglement entropy 
and has been proposed as a natural experimental way 
to measure many-body entanglement 57-59 . However this 
equivalence, which turns out to be a proportionality re- 
lation, has been only rigorously proven for free fermions 
in equilibrium, and verified in a Luttinger liquid with 
short-range interactions 59 . In particular for these two 
cases J-(x) ~ logx. Nevertheless, the equivalence breaks 
down for the dynamics of these two quantities after quan- 
tum quenches. For example, after an interaction quench 
in a short-range Luttinger liquid (a CFT), the entangle- 
ment entropy grows linearly with time until it saturates 
to a value that depends linearly on the size of the subsys- 
tem 24 , while the fluctuations T{x) grow logarithmically 
with time and saturate to logx. Notice that fluctuations 
grow in time after the quench in the same fashion as they 
grow, in equilibrium, with the subsystem size, i.e. log- 
aritmically. This is related to the conformal invariance 
of the theory and the fact that it posses a dynamical 
exponent z — l 60 . 

In the Coulomb LM in equilibrium J-(x) ~ log 1 ^ 2 (x), 
grows slower than logaritmically suggesting that the en- 
tanglement entropy scales in a different way for systems 
with long-range interactions. Instead, after a sudden 
change in the range of the interaction, we find that the 
particle-number fluctuations behave as: 



Fix a-J^W' fi»*< 
J-W-S ^log(f), fort» 5 ^ TT , 



(55) 



where Vpt 3> d. They grow logarithmically with time 
and saturate to a value that depends logarithmically on 
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the subsystem size, i.e. they behave as in a quench in the 
short-range LM. The only manifestation of the Coulomb 
potential is reflected in the timescale at which T satu- 
rates. 



C. Turning off the long-range interactions 

Next we briefly consider the opposite situation to the 
one analyzed above, namely the case where initially the 
system is in the ground state of the LM Hamiltonian Eq. 
(1) with Coulomb interactions, but the time evolution is 
conducted by a Hamiltonian with contact interactions, 
characterized by parameters Kf and Vf. Using the gen- 
eral form Eq. (B2) adapted to this case, we find the 
following expression for the one-particle density matrix 
in the steady state at leading order: 



C^ r (x » d, t = oo) ~ e 3 ^ 



(K 2 + l)log 3 / 2 (x/d) 



(56) 



This result is completely equivalent to one in which 
the Coulomb interaction strength g is suddenly changed 
(suitably modifying the prefactors in the exponent) , since 
the correction factor in the long-times regime has the 
same form as the initial correlations. Also in the steady 
state, we find 



C^(x>d,oo) 

Cg (x > <f,00) 



a -£V2fl 1/! W a (x/«l) 



(57) 
(58) 



We observe that the three correlation functions have a 
similar asymptotic behavior for large distances and that 
they differ from the correlations in the ground state. In 
particular, we see that, remarkably, Cg does not contain 
the parameters of if/; it only depends on the initial state. 

Regarding the dynamics at short times, since the 
Hamiltonian driving the time evolution has a linear QP 
spectrum, the linear light-cone effect holds exactly. In 
fact, for times d <C %Vft <C x correlations look like those 
of the initial state up to a time dependent prefactor: 



C^ r (x,t) 



C£(x)e «r 



/-""(*)>' 



/2^lo e 3 ' 2 (2v f t/d) 

) 

r 2 V2 1/2 i 



(59) 



\x) e - 2m ^/^9 1/2 i°s 3/2 (2-/tM (60) 



Ci i),n (x)e^5 



g 1 ' 2 \og i ' 2 (2vft/d) 



(61) 



The initial condition is given by the equal-times correla- 
tion functions 33 , 



Cf' m {x) 
Cf"(z) 



x e 

-2m 2 



■ \og^ 2 (x/d) 



/2^1og 3 / 2 (x/d) 



(62) 

(63) 
(64) 



We notice that, as pointed out in the discussion on gen- 
eral quenches in Appendix B, if we completely turn off 



the interactions (Kf = 1) the one-particle density matrix 
has no time-dependence. 

Finally, the charge fluctuations follow the behavior: 



ff l/2 X 2 log 3/2 (f); for0> 



(65) 



2v f 



where Vft 3> d. Notice that the fluctuations grow in 
time in a complete different fashion (faster) as they grow 
in a quench starting from a short-range interacting ini- 
tial state (logarithmically). The saturation value is also 
different. This reflects the strong influence of the ini- 
tial condition on the subsequent dynamics of the charge 
fluctuations. 



IV. EXPERIMENTAL REALIZATION 

The experimental study of the out-of-equilibrium dy- 
namics of isolated quantum systems has become a reality 
thanks to the advent of systems of ultracold atomic gases 
loaded on optical lattices. This systems offer the possibil- 
ity of monitoring the dynamics of a quantum system that 
maintain coherence for long times compared to the typi- 
cal experiment durations as well as a high degree of tun- 
ability of the geometric configuration and the statistics of 
the particles (fermions, bosons or mixtures). Neverthe- 
less in the case of neutral atomic species, the interactions 
between the particles, although having a strength that is 
easily modified through Feschbasch resonances, are only 
of short-range (at most, they are of dipolar type, and de- 
cay as 1/x 3 ). In this section we propose ultracold ions 
trapped in electromagnetic cavities as a system suitable 
for the study of non-equilibrium dynamics of one dimen- 
sional systems of particles coupled by the Coulomb in- 
teraction. These systems have been intensely studied in 
the context of experimental quantum computation, since 
they satisfy the basic criteria for the implementation of a 
quantum computer (basically, they have long coherence 
times, the universal operations needed to implement any 
quantum algorithm can be realized using lasers tightly 
focused on the trapped ions, and their internal state can 
be accurately measured) 32 ' 61-64 . 

Ions can be cooled down to their ground state using 
laser cooling 65 and loaded on linear (ID) traps, com- 
monly called linear Paul traps. In the case of quan- 
tum computing experiments, the most popular setups use 
the ions |Be + (boson) or |oCa + (fermion) 32 , but a large 
variety of ions can be used. In these devices, a radio- 
frequency potential is applied to two electrodes which 
are parallel to the axis of the trap. These electrodes 
create an oscillating two-dimensional quadrupole poten- 
tial that is translational invariant along the trap axis. If 
the frequency of the radio-frequency field is sufficiently 
large, the ions experience an effective restoring force to 
the center axis. Additionally, static electric fields confine 
the ions along the trap axis. If the confinement perpen- 
dicular to the trap axis (radial direction) is much larger 
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t=0 



t>0 



FIG. 7: (Color Online) Scheme of the experimental realiza- 
tion. Initially the ions (blue) form dipoles with its image 
charges (dashed circles) placed symmetrically with respect to 
the gate (thick solid line). The dipoles interact with each 
other through the dipole interaction, decaying as x~ :i . Af- 
ter suddenly removing the gate the ions interact via the bare 
Coulomb repulsion. 



than the confinement along the trap axis, cold ions form a 
linear crystal 66 , in which the spacings are determined by 
a balance between the horizontal (axial) confining fields 
and mutual Coulomb repulsion. 

In the quantum computation setups 69 the measure- 
ments aim to the determination of the internal state of 
the qubit encoded by each ion. However, nothing pre- 
vents us from using the battery of techniques developed 
in the field of ultracold neutral atoms, such as the time- 
of-flight measurements, through which observables such 
as the momentum distribution and other few-point cor- 
relations are accessible on real time. 

The particular quench sequence proposed in this paper 
can be realized loading the ions on the linear trap in the 
presence of a gate parallel to the trap axis (see Fig. (7)). 
The interactions between the ions would be screened by 
the image charges generated by the gate, symmetrically 
placed with respect to it, leaving us with an effective 
dipole interaction between the atoms, in which case the 
system can be described as a short-range interacting Lut- 
tinger liquid. The sudden quench can be realized by 
rapidly removing the gate, allowing the ions to interact 
via the bare Coulomb repulsion. To reproduce exactly 
the sudden quench limit the gate should be removed in a 
time much shorter than any time scale of the system, or, 
equivalently, faster than any velocity scale. For a system 
of cold trapped loCa + ions, the Fermi velocity can be 
estimated from the typical densities 66 , and has an order 
of magnitude of lms -1 , but this can be reduced using 
heavier ions. However, based on the analysis made in 
Ref. 67, we can expect the sudden quench results to be 
valid for more general quench protocols, at least in some 
specific space-time regions. 



V. SUMMARY 



and vice versa. We have found that the stationary state 
correlations do not obey the typical power law decays, 
even in the case of turning off the Coulomb potential. 
This behavior can be obtained from a suitable defined 
GGE, which turns out to describe both static and dy- 
namic correlations. Moreover, the Coulomb potential 
modifies qualitatively and quantitatively some aspects of 
the dynamics while leaves other intact with respect to 
the case of a short-range interaction quench. In particu- 
lar, it substantially changes the light-cone effect picture, 
by introducing non linearities that mimic an acceleration 
of the excitations, and thus opening the possibility of 
a faster quantum information propagation. This effect 
can be understood based on the special characteristics 
of the plasmonic QP spectrum in ID, in particular the 
fact that the faster modes are the less energetic ones. Fi- 
nally, we have proposed that linear chains of ultracold 
trapped ions would provide a suitable setup for studying 
the effects described in this work. 



Acknowledgments 

We thank Miguel A. Cazalilla for enlightening discus- 
sions. This work was partially supported by CONICET 
(PIP 0662), ANPCyT (PICT 2010-1907) and UNLP 
(PID X497), Argentina. 



Appendix A: Formal solution 



In this Appendix we provide the solution to a gen- 
eral quench in the LM. We start by defining {a g ,al„} 



({< 



g, ' 



.}) as the diagonal basis of Hi (Hf). Since we 



shall work in the Heisenberg picture, the solution to the 
problem of the system's evolution can be expressed, at 
least formally, giving the time dependence of the opera- 
tors {a 9 ,al g }. We will use the explicit form of the Bo- 
goliubov transformations relating these operators to the 
free basis {b q , &!„} defined in section II: 



= M {l) 



(Al) 



and 



a_ 



(A2) 



We have presented a detailed study of the behavior with 
of the correlations following a special quench protocol in 

the LM, namely, to suddenly change the range of the in- \a( 1 ) — f cosn (^(<7)) smn (^(<7)) \ /A<n 



teractions between the fermions from short to long-range * I shxh(0/(g>)) cosh(#;(g)) 
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, where I — i, f. The transformation relating the initial 
and final basis is thus: 




_ { cosh(9 f (q)-9 l (q)) smh(0f(q)~9 l (q))\ ( a l \ 
\ wMfiM - 6i(q)) cosh(9 f ( q ) - 9M)) JW-J' 

(A4) 

Since the {a q , operators have trivial time depen- 
dence, i.e. a q (t) — a q exp[—ief(q)t], we arrive at the 
solution: 

a q (t) = f(q,t)a q +g*(q,t)at q (A5) 

where 

f(q, t) = cos(e/ (q)qt) ~ i sin{e f (q)t) cosh[2(0,(<f) ~ 9 % (<?))], 

(A6) 

g(q,t) =ism(e f (q)qt)smh[2(6f(q) - 0* (<?))]. (A7) 

Appendix B: General form of the correlation 
functions after an arbitrary interaction quench 

We provide here a general formula for all the corre- 
lation functions studied in this work in the case of an 
arbitrary quench. Defining as before: 

C 6 (x, t) ee (e m f t 6(x)6(0)e- iH f t ), (Bl) 

we find that 

C 6 (x,t) = C { ?\x)exp[-$ 6 (x,t)}, (B2) 

with 



* 6 (x,t) = -2^ e — l 6 {q)smh[2{6 f {q) - 9^))] 

q>0 q 

x (1 — cos(qx)) sin 2 (ef(q)t). 

(B3) 

The function ^^{x^) defines the correction introduced 
by the quench, (x) is the correlation function at t = 
0, i.e. the equilibrium equal-times correlation function 
for a system described by Hi, and 

! sinh(29 f (q)) if = ip r (x) 

2m 2 exp[-29 f (q)} if 6 = e 2lm ^ (B4) 
-^exp[26 f (q)\ if 6 = e in8( - x \ 

is a function depending only on the operator considered. 
It is worth noticing that if Jfi(q) = the correlation 
function has no temporal dependence, retaining the same 
spatial dependence as in the initial state. For example 
the one-particle density matrix and accordingly the mo- 
mentum distribution function are unaffected when the 
interactions are suddenly turned-off, since in that case 
sinh(26» / (<j)) = 56 . 

We shall illustrate this result following in some detail 
the calculation of C™(x,t). The case of Cg(x,t) is com- 
pletely analogous and the case of the Green's function 
has been already analyzed in detail in Ref. 23. The mode 
decomposition of the 4>{x) field is 




(B5) 

and using Eq. (A5) to introduce the time dependence, 
we find that 



({<P(x,t)-q>(0,t)} 2 ) =Y, £ e^ 2 ^ (9) (l-cos(^)){l+2sinh 2 (^<q)-^( 9 ))+sinh[2(0 / ( 9 )-^((7))]cos(2 e/ (g)t)}. (B6) 

If we separate the sum isolating a time independent part and a time dependent one that vanishes at t = we find the 
following factorization for C2 l (x,t): 

C?(x, t) = cf' m {x) exp{-4m 2 £ e — e -™^)(l _ cos ( qx )) s inh[2(^(g) - 9 t {q))] sin 2 (e g t)}, (B7) 

I 

where C^' m (x) — exp{— 2m 2 e^ 29i ^ J2 q >o ^TT~(^ ~ ground state of Hi, verifying the general result. 
cos(qx))} is exactly the equilibrium correlator in the 
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